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pH . A novel viewpoint, i.e., adaptive synchronization, is proposed 

to explore collective rhythm observed in many complex, self- 
organizing systems. We show that a simple adaptive coupling is 
able to tip arrays of oscillators towards collective synchronization. 
Two arrays of simple electrically coupled Hindmarsh-Rose chaotic 
■^j- ! neurons are used to illustrate cooperative dynamics of neural 

activity like the central pattern generators, which supplies a new 
t^. ; idea for biological experiments and numerical simulations. The 

results indicate that such small-world adaptive coupling may be a 
q ■ universal essence of the collective rhythm observed in the natural 

world. 
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Today one of the main unsolved problems in science is how to apperceive and study 
complex, self-organizing systems. A puzzling characteristic in these systems is the 
spontaneous collective rhythm, i.e., so-called collective synchronization. This re- 
markable phenomenon has been observed extensively in the natural world, ranging 
from inorganic systems to organic systems, e.g., Christiaan Huygens' two synchro- 
nization clocks, wobbly bridges, the oscillating uniformly Josephson junctions, the 
Bose-Einstein condensate (BEC), emerging coherence in chemical oscillators, fireflies 
synchronizing spontaneously their flashes, self-synchronization of the cardiac pace- 
maker cells, the cooperative pattern in network of neurons, animals' gaits, groups 
of women with the mutually synchronized menstrual cycles, an audience clapping 
in sync, and coherent moving states recently observed in highway traffic, etc^. In 
the collective rhythm, each individual itself executes two different behaviors, i.e., 
periodic oscillating and aperiodic (i.e., chaotic) oscillating. The collective rhythm 
of periodic oscillators emerges in many fields, while the collective chaotic synchro- 
nization is observed in the network of neurons and two weakly coupled BECs, etc. 
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As is well known, there exist infinitely many unstable periodic orbits embedded in a 
chaotic attractor, meanwhile due to the extreme sensitivity of initial values chaotic 
systems intrinsically defy synchronization. Therefore from this sense the collective 
synchronization of chaotic oscillators seems to be more difficult and significant than 
that of periodic oscillators. 

Although collective synchronization has been observed in various concrete ex- 
periments, some key theoretical problems remain open for scientists, especially for 
mathematicians. Polymath Norbert Wiener ever tried to develop a mathematical 
model of collection of oscillators in the late 1950s but little fruit was obtained' 2 ^. 
The theoretical breakthrough came from A.T. Winfree's pioneering work in 1966' 3 ', 
where he developed a mathematical model to study large populations of periodic os- 
cillators, and obtained some important results although stymied by the difficulty of 
solving the mathematical model. A crucial breakthrough came in 1975 when Y. Cu- 
ramoto refined and simplified the mathematical framework developed by Winfree^, 
where an exactly solvable mean-field model of coupled oscillators was proposed. 
Since these pioneering studies lots of theoretical works have been accomplished ' 5 l 
Note that only the case of coupled periodic oscillators, i.e., the collective phase 
synchronization, was investigated in the literature. Through the large amount of 
effort in idealized mathematical models there has been a common acknowledgement 
for the complex problem on collective rhythm: individual oscillators are coupled 
together, and there exists a critical value for the intensity of mutual coupling, below 
which anarchy prevails and above which coherence with the collective rhythm reigns. 
However, in this viewpoint two points deserve to investigate further. One is how 
the individual oscillators are coupled together, and the other is how complex, self- 
organizing systems produce the threshold of coupling strength to reach collective 
synchronization. In the previous theoretical work, the weak couplings were almost 
global, e.g., mean-field coupling, where each oscillator interacts with the rest of os- 
cillators. However from the perspective of biology (especially neurobiology) such 
all-to-all coupling is too complicated to benefit to realize promptly the correspond- 
ing function. As for the problem on the critical intensity of coupling, it is still far 
away from scientists. 

In my opinion, the architecture of interaction of weakly coupling in the collec- 
tion of oscillators should be as simple as possible, and the critical coupling strength 
should be reached self-adaptively due to the spontaneousness of collective rhythm. 
This speculation is inspired by the provocative words "self-adaption creates com- 
plexity" which is a basic idea of complex adaptive systems (CAS) theory, meanwhile 
this idea is motivated directly by considering the collective chaos synchronization. 
In network of chaotic oscillators, which emerges widely in the field of neurobiology, 
the mathematical model must be exactly unsolvable, but as is introduced above 
almost all the previously obtained results on collective synchronization of coupled 
oscillators are based on the solvability of the mathematical model. In this paper, 
a simple mathematical framework is proposed to confirm this speculation, where a 
simple adaptive coupling is used to produce the collective synchronization of arrays 
of oscillators. 



Our analysis will be limited to an array of oscillators that are all strictly identical, 
which ignores the diversity. Suppose that the dynamical behavior of each oscillator 
is governed by an n-dimensional differential equation 

x = f(x),xeK n , (1) 

where f(x) is a differentiate nonlinear vector function. We firstly consider cycle- 
type coupling in a ring of N oscillators (note that only the case of unidirectional 
coupling is considered here, but the extension of bidirectional coupling is very nat- 
ural), 

Xi = f(xi) - e M _i(xj - Xi-x), i = 1,2, •• -,7V, (2) 

where £j = (x^i, X^2, • • • , #i,n) £ ^0 = x Ni — Q,i-1,2, • - ' j e i,i-l,n), 

and €i t i-i{Xi — Xi-i) = (ej,j-i,i(xj,i — £1,1-1,2(^,2 — #1-1,2), • ■ • , £i,i-i,n{ x i,n — 

^i-i.n))- This model indicates that the connection topology is cycle-type, and each 
individual oscillator is coupled only to its nearest neighbor. Here the coupling 
intensity e^-i varies adaptively according to the following update law (this is very 
crucial) 

= li,i-l,j{ x i,j ~ x i-l,j) , i = 1, 2, ■ • ■ , N, j — 1, 2, • ■ • , 71, (3) 

where 7i,i-ij are arbitrary positive constants. Next we consider another simple 
configuration of network, i.e., star- type coupling, 

x 1 = f(x 1 ), ±i = f(xi) - €i tl (xi - xi), i = 2,3, • ■ -,N, (4) 

where all notations are as those above, and the coupling strength e^i varies according 
to the update law 

ei,ij = li,i,j{ x i,j ~ x i,j) 2 , i = 2,---,N, j = 1,2, ■••,n. (5) 

Such star-type coupling admits the characteristics of small-world networks, where 
the first oscillator, x±, represents a hub. 

Using the well-known Lasalle's invariance principle in mathematics, one may 
prove that for arbitrary bounded solutions of system (2)-(3) (or (4)-(5)) as t — > oo 
x i(t) ~ > %i{t),i = 2, • • • , N, and the coupling strength will converges to a constant 
dependent on the initial values. The idea of proof is similar to that in [6]. The results 
indicate that the above two coupling schemes can produce the collective synchro- 
nization of oscillators, and such collective synchronization is nonlinear stable and 
robust against the effect of noise. Moreover the variable coupling strength reaches 
self-adaptively a value which just corresponds to the threshold for arising collective 
rhythm, meanwhile the small weak coupling may be obtained by adjusting the dis- 
sipation parameter "fij-ij. In addition, note that in the present coupling scheme 
the mutual interaction will vanish once the collective synchronization is achieved, 
which implies that the collective rhythm doesn't change dynamical behavior of each 
individual. In particular, when the oscillator (1) is chaotic, thanks to the global 
attraction and nonhyperbolicity of chaotic attractor these coupling schemes will be 



more simple and the global collective synchronization will be realized more easily, 
see the following examples on cooperative dynamics in neurobiology. 

It has been observed in neurobiological experiments and in numerical simula- 
tions that individual neurons show chaotic spiking-bursting, while ensemble of such 
irregularly bursting neurons can produce cooperatively coherent, rhythmical burst- 
ing, i.e., synchronized pattern of neural activity' 7 '. Such cooperative property plays 
a crucial role in neural activity, for example the central pattern generators (CPG) 
controlling the rhythmic motor behavior of animals. An important question is how 
neural assemblies produce and control regular rhythm. In the literature, some global 
coupling schemes like mean-field coupling have been used to explore the interesting 
sync pattern l 8 '. Here we will attempt to apply the proposed adaptive synchroniza- 
tion to reveal this complex, self-organizing phenomenon. To do it, we choose the 
famous Hindmarsh-Rose neuron model t 9 ', a third-order ordinary differential equation 
with one slow variable modelling chaotic spiking-bursting of neuron, as dynamical 
equations of individual neurons. Firstly we consider a network of neurons as the 
following cycle- type electrically coupling scheme, 

Xi = Yi + SXf — Xf — Zi + I — e i) j_i(Xj — Xi_i), , . 

Yi = l- 5X? - Y, Z t = -rZ t + 4r(X t + 1.6), % = 1, 2, • • ■ , N, 1 ' 

with 

— li,i-l(Xi — Xi_i) , (7) 

where X = X N , r = 0.0012, and the external input I = 3.281. Each neuron is 
characterized by three time-dependent variables: the membrane potential Xi, the 
recovery variable Yi and the slow adaptation current Zj. e^j-i represents the self- 
adaptive strength of the electric coupling between the zth neuron and the (i — l)th 
neuron. Similarly, the star-type coupling network of neurons with the first neuron 
being the hub is given as 

Xi = Yi + SXf — Xf — Zi + / — 6i : i(Xi — X\), , s 

Y = l- 5Xf - Y, Z t = -rZ t + 4r(X, + 1.6), i = 1, 2, ■ • ■ , N, U 

with en = 0, and otherwise 

ei,i = 7i,i(*i-*i) 2 > i = 2,3,---,N. (9) 

For the sake of simplicity, we set the parameters 7^ = 0.1 and = 5. The 
adaptive collective synchronization of the two networks of chaotic neurons is shown 
numerically in Figs. 1-4. Comparing the above two coupling schemes (see Fig. 4 
legends), we find that the adaptive synchronization in the star- type network with 
the characteristics of small-world is achieved more easily, which confirms again the 
strong synchronizability of small-world networks' 10 '. In conclusion, such adaptive 
collective synchronization is more simple than those based on the global coupling 
schemes in literature, which supplies a new idea for biological experiments and 
numerical simulations. Moreover, the result supports the well-known neural theory 
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of high-level brain function proposed by D. Hebb in 1949: Learning takes place by 
the adaptive adjustment of the conductances of the connections between neurons^ 11 '. 

As the above statements, besides the self-adaption of individual neurons a key 
reason why such simple synaptic coupling schemes (i.e., only coupling one variable 
representing the membrane potential) can produce so interesting collective synchro- 
nization of neural activity is the chaotic characteristic of neural dynamics. This 
result throws highly light to the well-known viewpoint that chaos is a necessary 
ingredient in life. Obviously the present adaptive collective synchronization can be 
generalized to the other networks, such as the multi-star-type networks and the 
higher-dimensional lattice networks. In particular, such adaptive coupling can be 
extended to the complex networks like small- world networks and scale-free networks 
I 10 ' 12 !. Combining the characteristics of small-world networks, such as quick signal- 
propagation and strong synchronizability, etc., we speculate that such small- world 
adaptive coupling may be a universal essence of collective rhythm observed in the 
natural world. Therefore we hope that this work will inspire further studies of 
complex, self-organizing systems. 
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Figure. 1. The collective synchronization emerging in the cycle- 
type neural network (6). The collective synchronization in the 
simply model with 5-coupled chaotic neurons is shown by calcu- 
lating synchronization error in system (6), where ex,i = Xi — 
X\, ey,i = Yi — Y\ and ez.i = Zi — Z\,i = 2,3,4,5 respec- 
tively denoting the synchronization errors in each variable con- 
verge to zero. Here the initial values of variables are set as 
(0.2, 3, 0.7, 0.1, 4, 0.8, 0.3, 2, 0.6, 0.1, 2, 0.7, 0.3, 4, 0.6). 
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Figure 2. The temporal evolution of the mutual synaptic cou- 
pling strength e^j-i, i — 1, 2, 3, 4, 5 in system (6)-(7). The mutual 
synaptic coupling strength between neurons tends eventually to 
a small constant in the course of arising the collective synchro- 
nization in Fig. 1. The evolution is governed by the adaptive 
law (7), where the dissipation parameter is uniformly chosen as 
7i,»-i = 0.1, i = 1,2,3,4,5, and the initial coupling strength is 
uniformly set as 0. 
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Figure 3. The collective synchronization emerging in the star- 
type network (8) with 5-coupled neurons. The collective syn- 
chronization in the star-type network of coupled Hindmarsh- 
Rose chaotic neurons is shown by numerically calculating syn- 
chronization error in system (8), where ex,i = Xi — X\, 
e Y ,i = Yi - Y\ and e z ,i = Z { - Zi,i = 2,3,4,5 re- 
spectively denoting the synchronization errors in each variable 
tend to zero. Here the initial values of variables are set as 
(0.2, 3, 0.7, 0.1, 4, 0.8, 0.3, 2, 0.6, 0.1, 2, 0.7, 0.3, 4, 0.6). 
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Figure 4. The temporal evolution of the mutual synaptic cou- 
pling strength ei t i,i — 2,3,4,5 in system (8)-(9). According to 
the update law (9) the mutual synaptic coupling intensity be- 
tween individual neurons and the hub neuron self-adaptively ap- 
proaches the threshold for producing the collective synchroniza- 
tion in Fig. 3, where the dissipation parameter is uniformly chosen 
as 7^1 = 0.1, i = 2, 3, 4, 5, and the initial coupling strength is uni- 
formly set as 0. Obviously here mean coupling strength is smaller 
than that in Fig. 2. 
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